Colloidal aggregates tested via nanoindentation and quasi-simultaneous 3D imaging 
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The mechanical properties of aggregated cohoids depend on the mutual interplay of inter-particle 
potentials, contact forces, aggregate structure and material properties of the bare particles. Owing 
to this variety of influences experimental results from macroscopic mechanical testings were mostly 
compared to time-consuming, microscopic simulations rather than to analytical theories. The aim of 
the present paper was to relate both macroscopic and microscopic mechanical data with each other 
and simple analytical models. We investigated dense amorphous aggregates made from monodisperse 
poly-methyl methacrylate (PMMA) particles (diameter: 1.6 /xm) via nanoindentation in combina- 
tion with confocal microscopy. The resulting macroscopic information were complemented by the 
three dimensional aggregate structure as well as the microscopic strain field. The strain field was in 
reasonable agreement with the predictions from analytical continuum theories. As a consequence the 
measured force-depth curves could be analyzed within a theoretical framework that has frequently 
been used for nanoindentation of atomic matter such as metals, ceramics and polymers. The ex- 
tracted values for hardness and effective elastic modulus represented average values characteristic 
of the aggregate. On base of these parameters we discuss the influence of the strength of particle 
bonds by introducing polystyrene (PS) between the particles. 



I. INTRODUCTION 

Suspensions of micro- and nanometer-sized particles 
(colloids) were intensively used to model atomic systems 
[iHS] . Increased lengthscales in case of colloids enable for 
optical observation with light microscopes and lead to 
slowed down kinetics compared to atomic systems. Cap- 
illary waves could be visualized in real space [T, the dy- 
namics of dislocation generation in weak colloidal crystals 
was studies on a single particle level [3 and experiments 
on colloidal glasses contributed to a better understand- 
ing of the glass transition O |6] . These studies benefited 
from universal theoretical descriptions applied to both, 
colloidal and atomic systems. In particular mode cou- 
pling theory, initially developed for the glass transition in 
molecular liquids turned out to be suitable for the mod- 
eling of dense colloidal suspensions [7| fS]. In many of 
these systems the mechanical behavior is of central inter- 
est, as this characteristic is experimentally conveniently 
accessible. 

Besides this fundamental scientific motivation there is 
a technological interest in the mechanical properties of 
colloidal suspensions and in particular aggregated col- 
loids. Colloids are frequently used as filler materials in 
various products of daily life. Depending on the appli- 
cation certain requirements related to morphology, abra- 
sion resistance, specific surface and tendency to agglom- 
eration must be met. These properties depend on the 
physio-chemical properties of the used materials as well 
as the technical control of the production process. Yet, in 
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many cases the particles are not produced as single pri- 
mary particles but rather as aggregates consisting of sev- 
eral primary particles. Additional mechanical treatments 
are required in order to obtain the primary particles or 
certain aggregate sizes. Generally, the influencing factors 
for breaking up the aggregates depend on their structural 
and material properties [9] as well as the grinding and dis- 
persing process [TOl |TT] . Especially the micro- mechanical 
properties such as the breakage energy as well as plastic 
and elastic deformation energy at the nanoscopic scale 
have to be considered [12] [13] . 

Ordered and disordered assemblies of colloids also at- 
tracted interest because of their specialized functionality 
and advantageous preparation in large scales and via self- 
assembly techniques. Examples are photonic colloidal 
crystals that show a great potential as optical elements 
[m [15] or surfaces that can be made super-hydrophobic 
by coating them with raspberry-like particles [16]. Be- 
sides their primary function again the mechanical prop- 
erties play a key role as these colloidal aggregates must 
render a certain mechanical stability for a convenient us- 
age in daily life. 

Such systems fundamentally differ from those men- 
tioned in the beginning. Aggregated colloids are not in a 
dynamic equilibrium but are kinetically arrested. Bonds 
between the constituent particles last long and hinder 
structural rearrangements of the particles. Hence, any 
model based on the diffusive motion of the particles can- 
not be applied. Nevertheless, a comparison to equivalent 
atomic systems such as glasses, metals or ceramics pro- 
vides other ways for a theoretical description. It is the 
aim of the present study to test the applicability of ana- 
lytical continuum theory for the mechanical modeling of 
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We chose nanoindentation as the mechanical test 
method. In a nanoindentation experiment a tip of defined 
geometry, spherical or triangular-pyramidal, is impressed 
into the specimen and retracted while the applied force 
is measured. Simple, yet quantitative analysis methods 
based on the continuum mechanical models of Boussinesq 
[17] and Hertz [TF were successfully applied to interpret 
nanoindentation data. Although the specimen was tested 
only very locally the extracted mechanical quantities like 
hardness and Young's modulus were in good agreement 
with those obtained from classical macroscopic testings 
fl9t i2Qj . This might be related to the fact that despite of 
indentation depths down to 0.25 /im sufficiently large en- 
sembles of atoms or molecules were tested to be treated 
as a continuum and obtain average material properties 

In case of colloidal aggregates validity and applicability 
of these analysis methods had not been tested before and 
a priory could not be assumed. Even if the indentation 
depth was increased to several micrometers the number of 
involved particles was effectively decreased by a factor of 
roughly 10^ compared to atomic systems due to the size 
difference between atoms and colloids. Hence, the ques- 
tion arises whether the amount of tested material, i.e. 
the number of tested colloids, is large enough to properly 
define average material properties of the aggregate. The 
analogy is complicated even further if obvious differences 
in the interaction mechanisms between the constituent 
elements in atomic and colloidal systems are taken into 
account. Like in the dispersed state particles interact via 
electrostatic, steric and van der Waals forces. However, 
as soon as the particle aggregate structural properties 
become relevant, too. Dense and in particular crystalline 
structures show partial elasticity [22 while fractal gels 
tend to be deformed plastically [23|. If associated re- 
organization effects come into play material and surface 
properties of the particles need to be considered. Rough 
particles lead to increased frictional forces ^24n26 that 
may favor rotational over sliding motions [27] . 

Although not directly detectable we accounted for 
these microscopic processes by complementing the bare 
force measurement with real-space imaging of the col- 
loidal films using confocal microscopy. With the help of 
particle localization ^28| and tracking algorithms [29^ the 
particulate structure of the aggregate as well as three 
dimensional (3D) displacements of all individual parti- 
cles were extracted during the whole indentation process. 
These displacements enabled for the calculation of the lo- 
cal strain tensor. Since this approach requires high 3D 
position resolution only few studies fSHl ^5 followed a 
similar approach. Yet, in these studies the results were 
successfully compared to continuum models, e.g. the nu- 
cleation dynamics of dislocation defects during the in- 
dentation of colloidal crystals [H [Sj • 

The paper is organized as follows. After a short repe- 
tition of the standard analysis procedure for nanoinden- 
tation data (section [ll| we introduce the materials and 
experimental methods (section |III|) used in this study. 



Then we will discuss the microscopic distribution of dis- 
placements and strains during the indentation process 
and compare it to theoretical predictions based on con- 
tinuum mechanics (section IV A). In sections IV B and 
IV C the statistical analysis of the nanoindentation data 
will be presented and average material properties will be 
defined. 



II. ANALYSIS OF NANOINDENTATION DATA 

A typical force-depth curve of a complete loading- 
unloading cycle in a nanoindentation experiment is 
shown in Fig.[l] Upon pressing into the sample the inden- 
ter deforms the sample both plastically and elastically. 
The applied force F as well as the contact area between 
indentor and sample, projected to the surface along the 
indentation axis, A increases steadily. Both quantities 
are related with each other in terms of the hardness H 
of the material. 



H = F/A 



(1) 



This relation for the hardness is mainly motivated from 
analogous macroscopic indentation tests after Vickers 
^33j or Brinell ^34^ During the unloading process the elas- 
tic deformation is recovered resulting in a force pushing 
the indenter out of the specimen. This elastic recovery 
defines the reduced Young's Modulus Er via 
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The contact stiffness S is defined as the slope of the re- 
traction curve Fret{h) at the maximum depth /imax- 



(3) 



The Young's Modulus E of the sample is obtained by 
eliminating the influence of elastic deformations of the 
tip {Ei, 
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l/Er - (1 - V^)/Ei 



(4) 



Here v and Vi denote the Poisson ratios of sample and 
indenter, respectively. In contrast to the hardness the 
Young's modulus is well defined and equations ([2| to 
(|3| are based on the analytical continuum theories af- 
ter Boussinesq [17 and Hertz [TF. Originally derived for 
the indentation of a spherical tip in an infinitesimally ex- 
tended half space, it was shown that equations (2| to (|3| 
also hold for other indentor geometries [35l I36fr More- 
over, Doerner and Nix proposed that S is determined 
from a linear approximation to Fret at /imax [37]. More 
recently, however, Oliver and Pharr pointed out that the 
steady reduction of the contact area A while retracting 
the tip corrupts such a linear approximation [38^ . Instead 
a power-law behavior 
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FIG. 1. Exemplary force-depth curve: See text for explana- 
tions of used variables. (After Oliver and Pharr |41J). 



is more accurate and the fitting parameters hf and m 
allow for a calculation of S. 

Except for the projected contact area A all quantities 
needed for the evaluation of hardness H and Young's 
Modulus Er can be directly determined from the force- 
depth curve. In practice two routes are commonly ap- 
plied to measure A, imaging of the residual indent via 
electron, light or surface probe microscopy or a math- 
ematical procedure based on further assumptions about 
the elastic recovery. However, both approaches show ma- 
jor drawbacks and the method of choice depends strongly 
on the experimental details. 

The mathematical determination of A was proposed by 
Oliver and Pharr [38 and assumes a linear dependency of 
the restoring force on the indentation depth. Therefore, 
the effective depth he of the tip in the sample is given by 

he = /i-max — e i^max/ S (6) 

with a geometrical factor e. Once he is known A can be 
calculated from geometrical considerations [19], e.g. for 
a Berkovich indenter like in the present study: 

A = 24.56 hi (7) 

Beside the normalized mechanical quantities H and 
E'eff the total, elastic and plastic deformation works can 
be used to characterize the mechanical properties of the 
specimen [39", "40 . They are defined by the integrals of 
the indentation Find(^) and retraction curves Fret(^) as 
follows: 

/"^max 

Wtot = / Find{h)dh (8) 
Jo 

Wei.= / F,,t{h)dh (9) 
Jo 

Wpla = Wtot - Wela (10) 



III. EXPERIMENTAL 
A. Sample Preparation 

Colloidal aggregates were prepared as thin films. Com- 
pared to isolated aggregates such extended films are ad- 
vantageous as the mechanical properties are not affected 
by finite size-effects and results from different indentation 
locations are well comparable. 

The base particles were made from Poly(methyl 
methacrylate) (PMMA) and synthesized via dispersion 
polymerization with heptane as solvent [42] . A comb-like 
graft copolymer with a backbone of methyl methacrylate 
and glycidyl methacrylate and teeth made from poly (12- 
hydroxystearic acid) [43] was copolymerized with the 
PMMA and preferentially allocated at the surface of the 
forming particle. The poly(12-hydroxystearic acid) teeth 
served as a steric stabilization and prevented aggregation 
of the particles. Morevover, Nile red was added to the re- 
action mixture and homogeneously incorporated into the 
particles. The particles had a mean diameter of 1.6 jam 
and a polydispersity of 5 % obtained from scanning elec- 
tron microscope images on the basis of 100 particles. 
In order to minimize differences in the particle proper- 
ties, all used particles originate from one single synthesis 
batch. 

For a first set of films these PMMA particles were 
dispersed in hexane (Sigma- Aldrich, 98 % purity) at a 
low volume fraction of 2%voi. A 10 /il large droplet was 
placed and dried on a preheated glass substrate with 
150 fim thickness [44 at an elevated temperature of 50 °C 
without any further environmental control. Due to the 
fast evaporation and presumably strong convection dur- 
ing the drying process the resulting structure turned out 
to be completely amorphous. 

For a second set of samples a mixture of cyclohexyl 
bromide (CHB) and cis-decahydronaphtalin (DEC) was 
used to disperse the PMMA particles at a volume frac- 
tion of 10%voi- Various amounts of polystyrene (PS, 
64kg/mol, M^/M^ 1.03) were dissolved in 
the suspension to incorporate a tunable amount of solid 
bridges between the particles in the dried film. The mix- 
ing ratio of CHB and DEC was set to 4:1 by weight for 
practical reasons: Pure CHB swells and partially dis- 
solves the particles while films casted from pure DEC 
were not stable and peeled off the glass substrate af- 
ter evaporation. Moreover, the film surface was very 
smooth which might be related to the approximate den- 
sity matching of the mixture and PMMA. In the dis- 
persed state the dissolved PS acted as a depletion agent 
and induced attractive forces between the particles [45] . 
As a result the particles aggregated to extended clusters 
[46] . Yet, these attractive forces were small compared to 
the hydrodynamic forces acting on the particles during 
drying at 50°C. In the late stage of the drying pro- 
cess the remaining amount of CHB and DEC allocated 
on the surface of the particles forming capillary bridges 
between them. The PS was strongly enriched in these re- 



4 




FIG. 2. Scanning electron microscope images of the surface 
and interior a fractured film made of PMMA particles with 
additional PS (1 wt% with respect to the amount of PMMA). 
The homogeneous distribution of PS throughout the whole 
film is visualized in the large image. The PS coverage of the 
PMMA particles and the solid state bridges are revealed if a 
particle is broken out of the film (see inset). 
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FIG. 3. Radial distribution functions (RDF) for the two types 
of amorphous colloidal films, without and with PS added to 
the PMMA particles. The distributions are rather similar 
showing the typical oscillation behavior. The normalization 
constant do = 1-6 /im was the same for both structures. Com- 
pared to the films with PS those without PS the first peak 
dropped in amplitude and slightly broadened. This indicates 
an overall smaller volume fraction of particles and local den- 
sity variations, respectively 08] • 



gions forming solid glassy bridges (glass transition tem- 
perature of PS: 95 °C) when the solvent was completely 
evaporated. These bridges as well as the homogeneous 
average distribution of PS through the whole sample can 
be seen in the scanning electron microscope images in 
Fig-ll 

Further insight into the three-dimensional film struc- 
ture wa s gain ed with the help of confocal microscopy (see 

Infiltration of 



TABLE I. Preparation methods and structural parameters of 
the colloidal films. The errors give a measure of the variations 
for different imaged regions of the sample. 



section 



IIIB for additional information), 
the film with a liquid that matches the refractive index 
of the PMMA significantly increased the image quality. 
We used so called laser liquids with a specified refrac- 
tive index which are commercially available (refractive 
index: n = 1.5780 ±0.0002, Cargille Laboratories, USA). 
The films were neither re-dispersed nor were the solid 
state bridges affected which was checked up to several 
months after infiltration. The optical resolution was good 
enough to individually resolve the constituent particles. 
The thickness of the films ranged from 25 jam to 40 fim 
with local height variations of less than 2 jam over 10 /im 
in lateral dimensions. Using particle localization algo- 
rithms the 3D coordinates of all particles in the imaged 
volume were extracted. From these coordinates the ra- 
dial distribution functions (RDF) [47 of both types of 
films were calculated (Fig. [3|. Both distributions are 
very similar showing a strongly damped oscillation that 
is typical for an amorphous structure. However, a de- 
crease and broadening of the first peak at d/do = 1 was 
observed when introducing PS to the films. These vari- 
ations are most likely caused by a decrease in the total 
volume fraction and by local variations in the particle 
volume fraction [48j . 

The exact details on the particle interactions in the 
films were not investigated. The following assumptions 
were sufficient for the following discussion of the results. 
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fast drying 
at 50 °C 


amorphous 


(59 zb 3) 



in 4:1 (vol) mixture 
of CHB and DEC 



For the samples with added PS the particles predomi- 
nantly interact with each other via solid bridges. These 
bridges have high elastic modules (~ GPa) but are brit- 
tle and thus break upon small deformations. For the 
other samples dried from the suspension in hexane short- 
ranged contact forces of the small stabilization layer on 
the particle dominate the particle interaction. This inter- 
action again defines a threshold value for bond rupture 
which yet is much smaller than that of the films with 
incorporated PS. 



B. Confocal Microscopy 

We used a home-made laser scanning confocal mi- 
croscopes (LSCM) in fluorescence mode which is de- 
scribed in detail in reference The scan rate was 
one 2D frame per second with a sampling resolution of 
0.1/im/pixel in x and y-directions and 0.25 /im/pixel in 
z-direction. The microscope was equipped with an oil im- 
mersion objective (Olympus, UPlanApo PH3, NA=1.35) 
resulting in a spatial resolution of 250 nm in lateral direc- 
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tions (x, y) and 550 nm along with the optical path (z) 
if the refractive index was matched as mentioned above. 
Further general information about LSCM [50l |5T] as well 
as the special case of fluorescent labeled particles [52 [53] 
can be found elsewhere. 

The determination of the particle coordinates from the 
3D microscope images was done using the algorithm by 
Crocker, Grier [28 and Weeks [6 . In the first step the 
image quality was enhanced by a spatial filtering using 
a Gaussian kernel. Afterwards the particles appeared as 
isolated local intensity maxima which could be located 
subsequently. An additional position refinement step in- 
creased the precision of the particle coordinates to a level 
of better than 0.03 /im. At this point the image analy- 
sis was finished and the colloidal film structure was well 
defined by a set of 3D particle coordinates. 



Nanoindentation 




particulate y continuum 
z 



FIG. 4. Indentation geometry: The spherical indenter tip is 
indented in positive z-direction with h denoting the depth of 
the tip inside the sample. Due to the cylindrical symmetry 
of the setup only the xz-plane is shown. We refer to the 
colloidal film either in the native particulate picture r ^^-^ or via 
continuous position vectors r. Next to Cartesian coordinates 
also cylindrical coordinates i^, Lp are introduced. 



A simultaneous microscopic deformation analysis and 
force measurement was not possible with the present 
setup. Therefore, two complementary experiments were 
performed, force-sensitive nanoindentation with a com- 
mercial nanoindentor and so called "live indentation" in 
combination with confocal microscopy. 



Force- sensitive nanoindentation 

Force-sensitive nanoindentation experiments were 
done with a Tribo-Indenter® TI 900 from Hysitron Inc. 
equipped with a Berkovich diamond tip. We chose a 
displacement controlled operation for which the tip was 
pressed into the sample at a constant indentation rate. 
Besides the depth of the tip inside the sample the nor- 
mal force on the tip was recorded for the loading and 
unloading cycle resulting in so called force-depth curves. 
The maximal indentation depth and indentation velocity 
were kept constant for all samples at 3 /im and 0.5 /im/s, 
respectively. In order to get an estimation of the stan- 
dard deviation of the mechanical properties each sample 
was tested on a 5 x 5 pattern of individual indentation 
spots that were separated by at least 30 /im. After inden- 
tation the films were infiltrated with Cargille laser liquid 
and subsequently imaged with the confocal microscope in 
order to determine the extension of the residual indent. 



Live indentation 

For the live indentation experiments an indenter tip 
was mounted to a piezo-translation stage (Actuator PXY 
200SG, Controller ENV40, piezosystems jena GmbH, 
Germany) which in turn was placed on the sample stage 
of the confocal microscope. The tip was pressed into 
the film in steps of (0.5 — 1) jam. After each step the 



whole structure was imaged. During imaging the inden- 
ter was not moved. With the help of tracking algorithms 
[29 the trajectories of the particles could be followed 
through all measured 3D images. Thusly we obtained 
the vectorial displacement for each individual particle in 
the film for the complete indentation process. Besides a 
Berkovich tip also a spherical tip made from silica (diam- 
eter: 25 /im) was used. The primary purpose of the live 
indentation experiments was a comparison of the defor- 
mation field with the analytical theory according to Hertz 
^18^ |54 and Huber [55 . As this task was accomplished 
much easier for cylindrically symmetric indentation ge- 
ometry we will focus on the results using the spherical tip. 
However, essentially the same observations were made for 
the Berkovich indenter as well. As the imaging and in- 
dentation occurs quasi-simultaneously the film needed to 
be infiltrated from the beginning. 



D. Microscopic deformation analysis 

In case of nanoindentation of atomic systems the ap- 
plied external forces are transmitted through the speci- 
men resulting in stress and deformation fields that are 
described as continuous fields. The relationship be- 
tween these quantities is given by material properties like 
Young's and shear modules using the continuum elastic- 
ity theory, see e.g. Hertz [18, 54 and Huber [55 . 
On the contrary in particulate systems stresses a^*-^^ are 
localized to the contact points between two adjacent par- 
ticles i and j (Fig. [4| and a priori need to be described 
discretely. Depending on the type and strength of the 
particle bonds these stresses result in a change of the 
bond vectors d'^'^^\h) relative to the bond vectors before 
indentation d (0 /im) 

Ad^'^\h) = d^'^\h) - d^'^\0 ^m) (11) 
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with h denoting the depth of the indenter tip inside the 
film (Fig. [4|. Yet, every particle is directly influenced 
by all of its Nn nearest neighbors and particle rearrange- 
ments include rotations as well as tensile and compressive 
deformations of the local particle arrangement. We use 
the local elastic strain tensor e*^*^ [21 |3T] to quantify these 
deformations. It relates the changes in the bond vectors 
with the relative positions of the particles via an affine 
transformation. Any deviation from this affine transfor- 
mation is caused by irreversible reorganizations of the 
particle structure and errors in the particle positions. So, 
e^*^ can be determined by minimizing the expression 



a) 



Ad {h) - (h) d (^^'^ (0 Mm) 



(12) 



In principle Nn contains only particles, that are in 
physical contact with the central particle. However, poly- 
dispersity of the particles complicates the decision which 
particles actually are in contact. Yet, we checked the re- 
sults presented below for several values of the maximal 
distance between two neighboring particles and did not 
find a significant influence. 

Although e^*^ quantifies the deformation of the local 
particle structure it does not contain any information 
about the exact reorganization mechanisms. Elastic de- 
formations of the particle bonds compete with more com- 
plex reorganizations like rolling and sliding if the particle 
bonds break. Moreover, deformations of the particles it- 
self must be taken into account if the bonds are strong 
or local rearrangements are hindered. 
We also calculated the vectorial displacement 

(/^) ^ ^{i) (/^) _ ^{i) (0 ^m) (13) 

for particle i and its absolute value D^'^^ = ID*^*^]. As 
D*^*^ is based on relative particle positions, increased rel- 
ative errors are expected compared to the pure particle 
coordinates. This is even more pronounced for e that 
depends on changes in the relative positions. Therefore, 
both quantities were averaged over all particles that are 
located in a cubic subvolume with a side length of 3 jam 
which corresponds typically to an averaging over 8 par- 
ticles. Moreover, we only discuss data at the maximal 
indentation depth of h = 3 /im. In doing so we define 
quasi-continuous magnitudes D(r) and D{r) as well as 
e(r ) and J^{r ) at the position r. In the following we 
mostly show these averaged quantities in order to com- 
pare our findings with prediction of a continuum model 
of the indentation process. 



IV. RESULTS AND DISCUSSION 

A. Displacement and strain field of the amorphous 
structure 

We start the discussion of the results with the live in- 
dentation of the amorphous structure that was not rein- 
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FIG. 5. a) Distribution of absolute particle displacements 
D(r) in the xz-plane in case of an amorphous structure at 
the deepest indent position of 3 /xm. b) The average abso- 
lute displacement shows an exponential dependency on the 
distance to the indenter tip. Due to the resolution limit of 
about 0.05 /im a determination of particle displacements is 
not possible for distances exceeding 20 /jui. The position vec- 
tor r and the cumulative displacement vector D(r ) are mainly 
parallel as depicted in the inset in terms of the distribution 
of the projection r • D(r)/(|r| |D(r)|). 



forced by additional PS. As listed in table |l] the volume 
density was close to the value of 63 % expected for a ran- 
dom close-packed structure of monodisperse spheres. 
Fig. pk) shows the absolute cumulative displacement field 
D{r) for a xz-cut through the film at an indentation 
depth of 3 jam. The distribution was roughly cylindri- 
cally symmetric to the z-axis as expected for a spherical 
indenter tip and an amorphous structure. In this aver- 
aged representation no major spatial heterogeneities were 
apparent and the total displacement decayed monotoni- 
cally with increasing distance to the indenter tip. Fur- 
ther averaging over all azimuth angles (p G [0°,360°] 
and intervals of 15 ° in polar angles revealed an expo- 
nential dependency [see Fig. ^jp)] with no significant de- 
pendency on 1^. The error bars represent the statistical 
standard deviation from the averaging process. Despite 
rather large values of the standard deviation in the order 
of the displacement itself the exponential fit is very good 
for the complete resolvable distance range. At a distance 
of about 20 fim the displacement dropped below the ex- 
perimental resolution limit of the particle displacements 
of about 0.05 fim ^ ^20.03 jam corresponding to uncer- 
tainties in the particle positions mentioned above. This 



low value for the resolution limit suggests that the large 
standard deviations have a physical origin: they reflect a 
rather broad underlying distribution of displacements. 
The overall uniformity of the displacement field was also 
observed for xy-cuts (not shown) at various depths. Dis- 
placement and particle position vectors were mainly par- 
allel to each other. This is shown in the inset of Fig. [sJd) 
in terms of the projection r • D(r)/(|r| |D(r)|) defined 
as the normalized scalar product of both vectors. The 
histogram has a maximum at 1 corresponding to fully 
parallel position and displacement vectors. 

More detailed information could be obtained from the 
strain tensor e. Because of azimuthal averaging only the 
positive X-axis is depicted in the graphs of e at an in- 
dentation depth of 3 /im in Fig. [6| Just like the absolute 
displacement the absolute values for the strain compo- 
nents decreased with increasing distance to the indenter 
tip. Besides, there was a distinct patterning of positive 
and negative strains. In case of the diagonal elements 
{§)aa = ^aa positivc and negative values represent dila- 
tion and compression of the local structure in direction of 
Of, respectively. Directly below the indenter the structure 
was strongly compressed in z-direction and dilated in x- 
and y-directions. The opposite was the case for polar an- 
gles i9 >55° closer to the sample surface. Here the local 
particle arrangement was compressed in x-direction and 
dilated in y- and z-direction. Dilation and compression 
largely canceled each other leading to an almost vanish- 
ing relative volume change 
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Shear deformations also took place. The corresponding 
strain components in the xz-plane of Fig. [6] were dom- 
inated by the xz-component. It was maximal at polar 
angles of about 55 ° for which the absolute values of 
Cxx and Czz were minimal. 

We compared the spatial distribution of the strain ten- 
sor elements with a continuum model by Hertz [18 and 
Huber [55^ J6 If the indented material was homogeneous 
and isotropic and the deformation was purely elastic, the 
stress tensor a could be calculated analytically and re- 
lated to the strain tensor via Hooke's law 
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FIG. 6. Spatial distribution of the strain tensor elements e^^ 
at a depth of 3 /im for the case of a spherical tip indented 
in an amorphous film: The data is averaged for all possible 
values of G [0 °, 360 °] and therefore only shown for positive 
x-coordinates. 



-0.5 -0.05+0.005 0.05 0.5 



o 



o 

Oh 



10 



20 



strain / 
volume 
change 





10 20 10 20 
X position / (|Lim) 



FIG. 7. Analytical calculation of the strain field based on an 
isotropic, elastic solid. Compare to Fig. [6] Since the compo- 
nent Cxy vanishes, its graph is replaced by emax- See text for 
further details. 
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with E and v denoting the Young's modulus and Poisson 
ratio of the film. 

The results of this calculation are shown in Fig. [7| The 
overall agreement with the experimental data was best 
for a Poisson ratio of 0.5. This value represents an in- 
compressible material that does not show any relative 
volume changes Av as it also was observed in the exper- 
iments. Moreover, the general agreement of calculations 
and experimental data in terms of the principle division 
into positive and negative values for the individual strain 
components {e)ij is clearly visible. Yet, the relative decay 
of the strain components with increasing distance from 



the indenter tip was much steeper for the experimental 
data. 

This discrepancy might be associated with the presence 
of the glass substrate which is not considered in the cal- 
culations. Finite size effects cannot be excluded but are 
probably less important as no infiuence of the actual sam- 
ple thickness was observed experimentally. Another pos- 
sible explanation is the presence of plastic deformations. 
As the bonds between the particles are rather weak for 
the films without PS, yielding of the film, i.e. irreversible 
particles rearrangement, must be considered. Yielding 
causes a release of the stress inside the film and thusly 
reduces the range of stress field compared to the cal- 
culations. Following Tresca's criterion for yielding such 
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FIG. 8. The transition from a weak dependency of the maxi- 
mal shear strain on the distance to the indenter tip to a power 
law behavior with a slope of approximately -3 is a hallmark 
of an elasto-plastic solid. The gray data points for the larger 
distances are masked by the resolution limit of the strain anal- 
ysis. The complete xz-view is shown in the inset with a color- 
scale identical to Fig. |6] 

plastic deformations will first take place at the niaxirQurQ 
of the principal stress difference [56]. For a spherical in- 
denter the maximum is located at a depth of 0.63 a below 
the indenter tip with a denoting the radius of the pro- 
jected contact area (Fig. [4| [54]. Indeed such a localized 
yielding had already been found for other colloidal sys- 
tems such as crystalline assemblies of soap bubbles [22] . 
In our case the yielding point was close to the indenter 
tip even for the largest indentation depth (0.63 a ^ 5 jiia) 
and masked by the limited spatial resolution after aver- 
aging. Moreover, if the yield stress was small, yielding 
could already take place at smaller indentation depths 
and therefore even closer to the tip where it could not be 
resolved. 

Hence, the colloidal film should rather be described as 
an elasto-plastic solid. In this model the strain in the 
specimen is divided into two regions [571 EH] • At small 
distances from the indenter r the local stress is large 
enough to exceed the yield stress of particle contacts. 
Plastic deformations take place by means of particle re- 
arrangements and the rate of stress increase with further 
indentation decreases. At larger distances beyond the so 
called elastic-plastic boundary the local stress drops be- 
low the critical yield stress and the deformation is elastic. 
This region is characterized by a fast dependency of 
the stress. Moreover, the distribution is merely indepen- 
dent of ^. Experimentally this stress characteristics is 
best mapped by the maximal strain difference Cmax as it 
is a measure of the maximal stress difference [56]. If the 
elements Cxy and Cyz are neglected with regard to exz [59] , 
emax is defined via 

^max = Vi^xx - 6^^)2/4 + 6^^ (16) 

The experimentally determined spatial distribution of 
Cmax as well as the distance profile (Fig. |8| are in rea- 
sonable agreement with the described model. 

The comparison with the continuum model indicates 



that the measured strain tensor has elastic and plastic 
contributions. In the microscopic picture these contri- 
butions can be identified by rearrangements in the local 
particle structure. If the applied forces are small particle 
bonds are deformed elastically but when stressed beyond 
the yield point frictional sliding of particles occurs [24]- 
[26] . These processes are largely irreversible. Compressive 
strains might also have a contribution from deformations 
of the particles itself. This deformation x is maximal 
for the particles that are in contact with the indenter 
tip. Under the assumption that a typical maximal load 
of 1.5 mN (Fig. [T]) is distributed over roughly 25 parti- 
cles, Hertzian contact theory [181 El] predicts a particle 
deformation of: 

/ gp2 \ 1/3 

Here = (1 /impart + ^/Rind)~^ denotes the effective 
particle radius defined by the radii impart and Rind of 
the particle and indenter, respectively. £^part ^2.5 GPa 
is the effective Young's Modulus of the PMMA and the 
indenter tip material defined in analogy of equation ([2|. 
However, the corresponding contact pressure 

1/2 

p=^S;,rt(;^) « 500MPa (18) 

by far exceeds the compressive strength of PMMA (^ 
80MPa). So the particles are plastically deformed which 
manifests in a flattening of the particle surface at the 
contact point and the calculated deformation of 140 nm 
is overestimated. Hence, these particle deformations were 
present but only had a minor contribution to the strain 
close to the indenter tip in Fig. |6] As the force on 
the particles decayed steeply when going deeper in the 
film only the upper most particle layers were affected by 
plastic deformations. The lower particles showed pre- 
dominantly elastic deformations that could not be dis- 
tinguished from structural reorganizations from the con- 
focal images alone. These considerations are also in ac- 
cordance with the observation from above, that the film 
is incompressible. 

B. Hardness and effective Young's modulus 

The possibility to model the microscopic deformation 
of a the film with continuum mechanical models suggest, 
that the force-depth curves can be analyzed in the frame- 
work described in section |llj However, at a closer look a 
few difficulties arise. 

At first the force-depth curves need to be corrected 
for the deformation of the thin glass substrate. Upon 
indentation the substrate deforms just like the sample 
and thusly reduces the actual indentation depth. The 
higher the indentation force the stronger the reduction 
is. Indentation on a bare glass substrate showed that 
the deformation is highly elastic and an effective spring 
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FIG. 9. The confocal image shows the colloidal film after 
indentation. The two neighboring indent appears as an dark 
area devoid of bright particles. The triangle indicates the 
shape of the indenter tip. 



constant of 20N/mm was obtained from the unloading 
part of the indentation process. The deformation of the 
substrate was calculated from the applied force according 
to Hook's law and subtracted from the raw indentation 
depth. 

Further problems arise from the imprecise starting 
point of the force-depth curves. Upon approaching to 
the sample surface the tip might indent directly on top 
of a particle or into the interstice between them. After 
the initial contact with the tip the upper most parti- 
cles rearranged to fit the shape of the indenter. As only 
few particles were involved this first consolidation process 
had a diffuse mechanical response that differed strongly 
from one indentation spot to the other. Large strains in 
the vicinity of the indenter tip were revealed in live inden- 
tation experiments and eventually led to the two small 
humps in Fig. [T] in the early stage of the indentation. 
Given the fact that the particle radius of 0.8 jam was not 
far off the total indentation depth /imax these onset un- 
certainties gave rise to large errors in he and even more 
pronounced in A according to equations (|6| and ^ . The 
problem of an unclear total indentation depth was even 
more intensified when adhesion came into play and single 
particles were pulled out of the film and adhered to the 
indenter tip. These adhered particles could cause a pre- 
mature start of the subsequent indentation process and 
with this again gave another source of errors in /imax- 

We bypassed these difficulties by choosing the second 
possibility to determine the contact area A by imaging 
the remnant indent. However, this method disregards the 
elastic recovery of the material upon unloading. Yet, the 
live indentation experiments showed that only particles 
with small displacements at larger distances to the tip 
return to their initial locations. On the contrary parti- 
cles in direct contact with the indenter were largely irre- 
versibly displaced and mapped the tip geometry. Fig. |9] 
shows a xy slice of a complete 3D confocal data set show- 
ing two neighboring indent positions. The triangular 
shape of the Berkovich indenter is clearly visible and a 
determination of A is feasible. 

The histograms of the extracted values for hardness H 
and effective modulus £^eff of the amorphous films with- 



out additional PS are shown in Fig. [TO] Although the 
distributions are rather broad, they can be well described 
by a single Gaussian function each. The mean values are 
given by F = 0.025±0.010 GPa and E = 0.50±0.10 GPa. 
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FIG. 10. Histograms 



The same is true for the films with additional PS: the dis- 
tributions are monomodal. However, already the small- 
est amount of PS led to a sharp increase in hardness 
and elastic modules by factors of 5 and 2.5, respectively 
(FigjTT]). Remarkably, additional PS did not cause any 
further enforcement of the film. 

Due to the lack of other approved mechanical test 
methods we can compare the observed mean values only 
with the corresponding literature values for bulk PMMA 
(i^iit = 1.9 GPa [60] and E^t = (1.8 - 3.1) GPa). The 
hardness of a material is defined as the force per contact 
area that is needed to deform the specimen irrespective 
of the type of deformation. Irreversible rearrangements 
of particles are not distinguished from elastic deforma- 
tion of particles and bonds. In colloidal films without 
additional PS rearrangements of the particles are pro- 
moted by a small bond strength leading to a relaxation 
of stresses. As a result the force required to indent the 
film is largely reduced compared to bulk PMMA. The 
hardness is decreased by almost two orders of magnitude 
compared to the literature values. The situation is differ- 
ent in case of the films with PS where strong solid bridges 
between the particles hinder rearrangements. We stress 
on the fact that the interaction between the particles is 
still mainly based on van der Waals forces. Hence, the 
strength of the particle bond is a mere effect of the in- 
crease contact area between adjacent particles mediated 
by the PS. As a result the hardness is increased but still 
far below the value for bulk PMMA. 

The Young's modulus E is di measure of elastic recov- 
ery in the very beginning of the retraction cycle. In this 
sense it is directly correlated to the elastic energy that is 
stored in the deformation of the contact bonds and the 
particles itself. If the particles do not rearrange or de- 
form plastically during the indentation a larger amount 
of elastic energy can be accumulated and recovered when 
the tip is retracted from the sample. As these rearrange- 
ments are partially hindered for the films with additional 
PS due to the increased bond strength, E is larger than 
for those without PS. However, further addition of PS 
does not affect the mechanical properties of the enforced 
films. This behavior can be understood, if one assumes 
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FIG. 11. Hardness and Young's modulus show the same be- 
havior for a variation of binding strength between the particles 
in the film. Upon a sharp increase for the smallest amount of 
added PS, both similarly approach a plateau. 



this relation metals, ceramics and glasses could be uni- 
versally superimposed on a single curve. The same was 
true for our colloidal films. Films with and without PS 
collapsed to a single line with a slope of (—0.35 ± 0.03) 
(Fig. 12). Even semi-crystalline colloidal films, which will 
be discussed in detail in further publications, follow the 
same relation. The values for the slope differs substan- 
tially from the corresponding value reported in literature 
of 0.2 [41 for atomic systems. This, however, might be 
caused by systematic errors in the values for the hard- 
ness and Young's modulus as possible sources for these 
uncertainties were discussed above. Only further exten- 
sive experiments or simulations can uncover the exact 
origin of the deviation. 
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FIG. 12. Irrespective of structure or binding strength the 
mechanical properties of all samples can be universally related 
with each other in a plot of the ratio H/ E versus the fraction 
of plastic deformation work Wpia / Wtot • 



that E of the films with PS is effectively determined by 
the elastic properties of PMMA and PS themselves. As 
the Young's modules of PS and PMMA at room temper- 
ature are in the same order of magnitude E of the film 
depends mainly on the total volume fraction of PMMA 
and PS together. This quantity in turn is supposed to be 
rather independent of the relative amount of PS since PS 
and PMMA distribute homogeneously during the prepa- 
ration of the films. This model is further supported be 
the fact, that E for the films with PS is very close to the 
literature value of bulk PMMA. 



C. Universality 

Further checks for the consistency of the used analysis 
procedures were found in the previous works [61] that 
related the mechanical quantities H and E to the elas- 
tic and plastic fractions of the deformation work. Scal- 
ing arguments predicted a linear correlation of H/E and 
Wpia/Wtot irrespective of the used materials ^39l|62]. Via 



V. CONCLUSIONS 



We presented a study of the mechanical properties 
of colloidal aggregates tested via nanoindentation. Al- 
though the sample was tested only very locally the in- 
duced deformation was in qualitative agreement with the 
predictions of a continuum theory as well as with inden- 
tation of atomic matter. As a result the theory of Oliver 
and Pharr was successfully applied to obtain average ma- 
terial properties. Unfortunately, there is a lack of other 
methods to verify the absolute values. Hardness and ef- 
fective elastic modulus showed a strong dependency on 
the strength of the particle bonds as it has strong im- 
pact on the dominant deformation processes during the 
indentation. 

In general nanoindentation proved to be a potent 
method for mechanical characterization of colloidal films. 
Besides a simple experimental procedure the measure- 
ment data can be processed with standard analysis meth- 
ods. Despite these benefits and the combination with 
quasi-simultaneous confocal microscopy no information 
about the microscopic forces between the particles and 
the exact reorganization processes could be gained. For 
this reason also the origin of the large statistical scatter 
in the average displacements of Fig. |5]b) or the mechan- 
ical quantities in Fig. 10 remains unclear. Whether a 



broad distribution of contact forces like in force chains 
can account for it needs to be clarified. 

We will complement our findings with discrete element 
method (DEM) calculations and further experimental 
studies in future publication. The DEM calculations can 
be used to simulate the microscopic deformation field as 
well as force-depth curves on the base of the mechanical 
properties of the constituent materials and thus is com- 
plementary to the presented experiments. On the other 
hand, in future experiments we will use hollow spheres as 
local force sensors [63] and anisotropic labeled particles 
for sensing individual particle rotations [49 in order to 
get a deeper insight into the reorganization processes. 
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